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1.  SUMMARY 


Subauroral  magnetic  field  lines  map  to  the  inner  magnetosphere,  including  the  ring  current 
(RC),  outer  radiation  belt  (RB),  and  plasmasphere  earthward  of  the  electron  plasma  sheet 
boundary  (the  auroral  boundary  in  the  ionosphere).  During  space  storms,  the  inner  edge  of  the 
RB  electrons  follows  the  perturbed  plasmasphere’ s  boundary  (the  plasmapause,  PP),  while  the 
plasma  sheet  (PS)  electron  and  ion  boundaries  significantly  separate.  Energetic  ions  form  the 
perturbed  ring  current.  In  turn,  channels  of  intense,  poleward  directed  electric  fields  develop  in 
the  subauroral  geospace  creating  jets  of  enhanced  westward  convection  termed  Sub-Auroral  Ion 
Drifts  (SAID)  and  Sub-Auroral  Polarization  Streams  (SAPS).  Observations  indicate  SAID/SAPS 
as  the  most  prominent  feature  of  the  perturbed  subauroral  geospace.  The  subauroral  region  is 
among  significant  contributors  to  space  weather  due  to  distinct  electromagnetic  and  plasma 
structures  and  turbulence  within  SAPS  [1,2].  These  result  in  precipitating  RC  protons  and  RB 
electrons,  high  TEC  gradients,  and  strong  scintillations  of  UHE  and  L-band  signals.  Besides  mid¬ 
latitude  scintillations,  the  SAPS-type  subauroral  processes  are  also  closely  tied  to  generation  of 
storm  enhanced  densities  and  creation  of  polar  cap  patches,  which  are  the  primary  driver  of 
scintillation  in  the  high-latitude  region.  The  development  of  subauroral  spatial-temporal 
variability  is  among  the  hottest  space  weather  topics. 

This  research  aims  at  attaining  the  required  level  of  understanding  of  the  SAID/SAPS  generation 
mechanisms.  The  objectives  are  to  characterize,  and  understand  basic  plasma  processes  that 
control  the  subauroral  geospace  dynamics/structure.  These  include:  (1)  the  interaction  of 
substorm- injected  plasma  with  the  plasmasphere  that  forms  a  plasmaspheric  boundary  layer  and 
creates  SAID/SAPS,  (2)  the  generation  of  SAID/SAPS-related  plasma  waves,  and  (3)  the 
development  of  a  numerical,  first-principle  hybrid  model  of  interchange  and  Kelvin-Helmholtz 
instabilities  in  the  midlatitude  ionosphere  which  includes  finite  Earmor  radius  effects.  The  results 
advance  the  understanding  of  the  physics  underlying  the  spatial-temporal  development  of  the 
perturbed  subauroral  geospace  and  pertinent  space  weather  effects. 

2.  INTRODUCTION 

During  last  three  decades  the  paradigm  of  voltage  and  current  generators  based  on  the  test 
(single)  particle  approach  dominated  the  subauroral  geospace  research.  The  former  places 
SAID/SAPS  between  the  inner  boundaries  of  >l-keV  ion  and  electron  convection  defined  by  the 
ExB  and  gradient  and  curvature  drifts.  However,  this  approach  violates  the  plasma  quasi¬ 
neutrality,  intrinsic  of  slow  processes.  The  current  generator  explores  closure  of  the  Region  2 
field-aligned  current  (EAC)  through  a  low-conductive  ionosphere,  resulting  in  locally  enhanced 
poleward  electric  field  Es.  The  generators'  concept  have  come  to  be  called  the  SAID/SAPS 
paradigm.  It  predicts  that  the  ring  current  and  SAID/SAPS  are  built  upon  a  slow,  a  few  hours, 
timescale.  The  corollary  to  the  above  is  that  SAID/SAPS  should  follow  the  asymmetric  RC 
development  and  be  enclosed  by  the  R2  and  R1  EACs.  However,  our  earlier  results  [3] 
demonstrate  that  this  paradigm  is  in  serious  error  and  cannot  explain  the  spatial-temporal 
variability  in  the  region  of  interest. 

One  of  the  most  striking  discrepancies  is  that  substorm  subauroral  events  appear  in  ~10  min  or 
less.  Eigure  1  exemplifies  such  events  from  the  DMSP  E14  satellite  that  appear  in  less  than  10 
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minutes  after  the  onset  of  the  substorm  expansion  phase  [4,  5].  From  top  to  bottom  are  the 
downcoming  electron  and  ion  differential  fluxes  in  the  range  of  0.03-30  keV,  horizontal  plasma 
drifts,  electron  and  ion  temperatures,  and  plasma  densities. 
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Figure  1,  DMSP  F14  observations  of  fast  time  (left)  SAPS  and  (right)  SAID  events 

As  the  fast  timescale  is  characteristic  of  propagation  of  substorm  injection  fronts,  we  focus  on 
the  observations  near  substorm  onsets  (see  sections  3.1  and  4.1). 

Ionospheric  scintillation,  a  key  component  of  the  Space  Weather,  disrupts  communication, 
navigation,  and  surveillance  systems  [6].  Scintillations  in  the  UHF  and  GPS  L-band  channels  are 
caused  by  small-scale  ionospheric  irregularities  with  transverse  scale  lengths  <100  m.  Strong 
UHF/GPS  scintillations  observed  within  the  SAPS  region  indicate  the  presence  of  such 
irregularities  [2].  As  the  state-of-the-art  fluid  models,  e.g.  AFRL/PBMOD  [7],  are  incapable  of 
predicting  small-scale  structures  less  than  1  km,  it  is  necessary  to  go  beyond  the  fluid 
approximation  and  include  ion  kinetic  effects.  Therefore,  we  develop  a  numerical,  physics-based 
hybrid  model  which  includes  finite  Larmor  radius  effects  and  thus  capable  of  describing  small- 
scale  irregularities  (see  sections  3.2  and  4.2). 

3.  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 
3.1  Selection  of  events 

The  Kyoto  WDC  quicklook  auroral  indices  AU/AL  and  AE  and  the  IMAGE  database  [8]  are 
analyzed  to  specify  the  development  of  substorms.  Then,  we  use  the  AERE  database  of  DMSP 
and  GRACE/CHAMP  and  the  NASA  Goddard  database  of  Cluster,  Polar,  and  CRRES  to  obtain 
key  ionosphere/thermosphere  and  inner  magnetosphere  data,  respectively.  The  main  initial  tusk 
is  to  identify  events  near  the  magnetic  equator  following  the  onsets  of  substorms  and 
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magnetically  conjugate  with  those  in  the  ionosphere.  After  such  events  are  identified,  they  are 
explored  to  establish  the  SAID/SAPS  salient  features. 

3.2  Tuning  a  hybrid  code  with  ion  kinetic  effects  to  ionosphere 

The  AFRL  high-resolution  code  "Flute"  [9]  was  originally  aimed  at  interchange  instabilities 
in  high-beta,  high  energy  plasmas.  To  apply  it  to  the  midlatitude  ionosphere,  a  system  of 
nonlinear  equations  is  derived  for  description  of  the  plasma  instabilities  resolving  spatial  scales 
down  to  '-'10  m,  as  needed  to  describe  scintillations  of  the  UHF  and  GPS  L-band  signals. 
Numerical  solutions  provide  spectra  of  developed  turbulence  for  various  combinations  of  the 
known  driving  forces,  including  velocity  shears,  neutral  winds,  and  penetrating  electric  fields. 
Validation  of  the  results  starts  comparing  with  the  large-scale  fluid  modeling  by  turning-off  the 
kinetic  effect.  The  ultimate  validation  of  the  model  should  be  achieved  by  comparing  the  high- 
resolution  spectra  with  those  revealed  from  the  high-resolution  in  situ  observations. 

4.  RESULTS  AND  DISCUSSION 

4.1  Observations  of  substorm  SAID/SAPS 

Overall,  more  than  40  substorm  events  in  both  hemispheres  with  the  time  span  after  the  substorm 
onset  ranging  from  10  min  to  2  hours  and  during  major  storms  have  been  explored  to  establish 
the  SAID/SAPS  salient  features  [10-13].  Figure  2  exemplifies  the  SAID  events  we  focused  on 
during  this  work.  Shown  are  the  SAID  channels  after  the  onset  of  a  substorm  [8]  near  L  =  6.8  and 
22.9  MLT  at  09:47:35  UT  on  18  March  2002.  Left  frame  displays  the  meridional  electric  field 
(not  to  scale)  from  the  Cluster  1  satellite  during  inbound  (the  southern  hemisphere)  and  outbound 
(northern)  passes  (Cl  at  an  altitude  /2-'25000  km)  in  the  meridional  plane  at  23.0  MLT,  with  the 
positions  of  Polar  (/2-'5500  km)  and  DMSP  F14  (/z-'850  km)  indicated.  Dipolar  magnetic  field 
lines  are  shown  for  reference.  The  thick  arrow  indicates  an  oncoming  substorm-injected  plasma 
jet  (plasmoid).  Right  frame:  (top)  Meridional  (outward)  electric  fields  in  mV/m;  (2"‘*row)  H+ 
differential  number  fluxes  (the  dashed  curve  shows  the  potential  energy  cAO  calculated  from  the 
measured  electric  field);  (3'^‘*)  Cl  1-keV  electron  counts  and  Polar/F14  electron  fluxes;  (bottom) 
the  plasma  densities.  Color  codes  are  in  logarithmic  scale. 

Apparently,  there  are  two  magnetically  conjugate  groups:  Cluster  -  F 14  in  the  southern 
hemisphere  and  Cluster  -  Polar  in  the  northern  hemisphere.  The  Cluster  and  Polar  electric  fields 
mapped  into  the  ionosphere  match  that  from  F14  to  within  ~10%.  The  hot  electron  population  at 
Polar  closely  resembles  that  mapped  from  DMSP,  and  the  1-keV  electron  flux  scales  to  Cluster's 
1-keV  electron  counts.  The  SAID's  outer  side  is  aligned  with  PP,  while  the  inner  edge  co-locates 
with  the  sharp  decrease  of  the  hot  ion  flux.  The  obvious  corollary  is  that  the  SAID  channel  is 
inside  the  plasmasphere  devoid  of  the  hot  PS  electrons.  Further,  the  test-particle  approach 
predicts  that  the  auroral  boundary  follows  drift  trajectories  of  (test)  keV  electrons  so  that  their 
maximum  energy  gradually  decreases  inward.  However,  this  disagrees  with  the  in  situ  data 
during  the  SAID  events.  Apparently,  the  hot  electron  boundary  is  coincident  with  the  density 
increase  on  the  plasmapause  and  represents  an  abrupt  dispersionless  cutoff. 
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Figure  2.  The  18  March  2002  Cluster/Polar/DMSP  SAID  event 


We  emphasize  that  this  type  of  the  substorm  auroral  boundary  remained  a  mystery  since  its 
discovery  [14],  Namely,  the  dispersionless  nature  of  the  hot  electron  cutoff  and  its  dependence 
on  the  cold  plasma  density  exclude  the  usual  (test  particle)  Alfven  layer  interpretation.  Nor  can 
the  test-particle  ’’nose”  distribution  of  westward-drifting  energetic  ions  fit  the  observed  hot  ion 
distribution  within  the  SAID  channel.  By  the  same  token,  the  generator  paradigm  cannot  explain 
the  size  of  the  SAID  channel  in  terms  of  the  Alfven  layer.  However,  these  features  follow 
directly  from  the  concept  of  a  turbulent  plasmaspheric  boundary  layer  (see  section  4.2). 

Figure  3  shows  the  related  magnetic  variations  (2"‘*  row)  and  currents  (3'^‘^  row)  for  the  SAID 
events  on  18  March  2002  and  8  April  2004  [cf ,  10,  Figures  6-8]. 


Figure  3,  SAID-related  current  systems 


Apparently,  the  downward  FACs  in  the  ionosphere  concentrate  near  the  SAID's  outer/poleward 
boundary.  In  the  magnetosphere,  the  downward  FACs  near  the  outer  edge  (^  =  0)  virtually  mimic 
ionospheric  counterparts.  However,  FAC’s  magnitude  at  Polar  is  close  to  that  mapped  from 
DMSP  but  twice  as  smaller  as  that  mapped  from  Cluster.  This  substantial  difference  and  the  lack 
of  the  ionospheric  counterpart  of  the  inner,  alternating  part  of  the  magnetospheric  current 
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indicate  partial  closure  of  the  downward  FAC  above  Polar  (>5500  km).  Note  also  that  no  density 
trough  is  present  in  the  ionosphere  (4*  row).  That  is,  eontrary  to  the  current  generator  paradigm, 
neither  the  field-aligned  currents  enelose  the  SAID  ehannel  nor  the  plasma  density  (conductivity) 
minimizes  in  the  SAID  region. 

The  overall  dataset  [3,5,  10-13]  shows  that  the  generator  paradigm  is  in  serious  error.  Therefore, 
we  put  forward  a  concept  of  the  SAID  channel  as  an  integral  part  of  a  turbulent  plasmaspherie 
boundary  layer  (TPBL)  formed  in  the  evening  seetor  where  the  plasmasphere  short-eireuits 
reeonnection-injected  hot  plasma  jets. 

4.2  Turbulent  Plasmaspherie  Boundary  Layer 

In  agreement  with  the  three-deeade  statisties  [e.g.,  15],  the  explored  fast  time  SAID  events 
eoncentrate  in  the  premidnight  MLT  sector,  characteristic  of  the  propagation  of  substorm 
injections.  We  consider  near-tail  substorm- injected  -  e)  plasma  jets  as  hot  plasma  jets  of  the 

cross-tail  width  a  moving  with  the  bulk  speed  Vo||x  aeross  the  magnetie  field  Bi  (x)  ||  z.  This 
system  of  reference  is  elose  to  GSM  near  the  magnetie  equator  and  midnight.  The  magnetic 
pressure  in  the  inner  magnetosphere  usually  greatly  exeeeds  the  ram  pressure  of  the  injeeted 
plasma  jets.  Accordingly,  infinite-width  jets  will  be  halted  near  the  source.  However,  jets  of 
limited  a  (hereafter,  plasmoids)  ean  propagate  across  magnetic  barriers  can  penetrate  magnetic 
barriers  owing  to  self-polarization  fields  formed  at  the  front.  For  that  their  density  Nh  must 
exeeed  Nmin~  0.3(Bi/400)^  cm“^,  where  Bi  is  in  nT  (see  details  and  referenees  in  [10]).  In  the 
events  of  interest,  a  typieal  upstream  hot  ion  density  Nh  ~  1-2  cm”^  amonts  to  ~(3-6)Nmin. 

However,  this  self-similar  motion  breaks  as  soon  as  the  surrounding  plasma  shorts  out  the 
polarization  field  at  the  front.  This  oecurs  at  the  plasmapause  when  the  cold  plasma  density 
exeeeds  the  critical  value  ~10  cm”^.  As  a  result,  the  hot  eleetrons  eome  to  a  halt,  thereby  forming 
the  dispersionless  plasma  sheet/auroral  boundary.  In  the  plasmasphere  devoid  of  the  hot 
electrons,  the  eold  plasma  maintains  quasi-neutrality.  This  is  a  short  eircuit  of  a  hot  plasmoid 
over  the  cold  plasma.  The  hot  ions  yet  move  inward  over  the  distance  ~(0.1-0.2)Re  until  being 
halted  by  the  emerging  SAID  eleetric  field. 

Therefore,  the  SAID  ehannel  can  be  pietured  as  a  turbulent  plasmaspherie  boundary  layer  with 
an  unstable  mixture  of  the  hot  ion  and  eold  plasma  populations.  The  eritieal  role  of  plasma 
turbulence  for  the  SAID  formation  is  to  provide  anomalous  resistivity  in  the  circuit  and  magnetic 
diffusion.  Save  the  short-circuiting  and  plasma  turbulence  effects,  quite  similar  approach  is  used 
for  understanding  impulsive  plasma  penetration  of  the  dayside  magnetopause. 

Let  explore  the  wave  aetivity  accompanied  the  events  in  question,  alike  the  SAPS  events  [e.g., 

1].  Figure  4  presents  electric  (mV  /m  Hz)  and  magnetic  (nT  /Hz)  wave  eomponents  from  Cluster 
eoincident  with  the  above  SAID  events  [cf ,  10,  Figure  4].  Horizontal  lines  on  the  frequency-time 
speetrograms  indieate  the  lower  hybrid  resonanee  (solid)  and  the  tenth,  fourth,  and  second 
harmonics  (dashed)  of  the  /7T-ion  gyro  frequency  ~5.4  Hz,  derived  from  the  observed  magnetic 
field.  It  is  seen  that  the  wave  activity  is  greatly  enhanced  inside  the  SAID  channel.  Intense 
broadband  waves  arise  in  the  entry  layer  of  the  thickness  '-90  km  near  the  outer  boundary  and 
extend  to  the  eenter.  The  waves  in  the  entry  layer  are  mainly  eleetrostatie  ion  cyclotron/Bemstein 
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(ICB)  and  lower  hybrid  (LH)  modes,  while  highly  oblique  eleetromagnetic  whistler  (W)  and  fast 
magnetosonie  (MS)  waves  appear  near  and  earthward  of  the  Es  peak. 


18-Mar-02 


8-Apr-04 


14151617181920  54  57  60  63 

minutes  after  1 0:00  UT  minutes  after  0600  LIT 


Figure  4,  SAID-associated  plasma  waves 


Figure  5  gives  more  details  of  the  18  March  2002  southern  SAID.  From  top  to  bottom  are  (a)  the 
C 1  outward  and  scaled  eastward  electric  components,  (b)  differential  omnidirectional  number 
fluxes  in  the  range  1-10  eV,  (c)  the  r.m.s.  amplitudes  of  the  electrostatic  waves  in  the  lower 
hybrid  (LH),  magnetosonie  (MS),  whistler  (W)  bands,  (d)  the  perpendicular  hot  ion  (electron) 
pressure  Pi^e,  density  Ni  and  ion  anisotropy  5T  =  Ti/T||-1,  (e  )  deviations  from  the  model  field  of 
the  observed,  AB  =  B-Bq  (the  thick  line),  and  the  'inner'  or  downstream,  AB/’^*  =  B/’'*  -  Bq 
(dashed),  magnetic  fields;  as  well  as  the  diamagnetic  correction  ABi+e*^'^ .  Here  ABd®’‘*  is  obtained 
by  a  polynomial  fit  of  the  magnetic  field  downstream  of  the  channel  (at  >10:18:10  UT)  and 
extrapolated  upstream. 

Upstream  of  the  channel,  we  have  taken  Pi+e=Pi+Pe  ~  (3/2)Pi  (the  hot  ion  pressure  from  Cluster) 
by  extrapolating  the  Polar/HYDRA  data.  The  excellent  matching  of  ABi+e‘^“^  with  the  observed 
AS  upstream  of  the  channel  (at  <10:15:10  UT)  indicates  the  permeation  of  the  incoming  hot 
plasma  by  the  downstream  field  (due  to  anomalous  magnetic  diffusion).  This  agrees  with  the 
value  of  the  effective  collision  frequency  due  to  wave-particle  interactions  calculated  with  the 
observed  turbulent  fields  [10].  However,  ABi+e‘^‘^  significantly  deviates  from  AS  near  the  outer 
boundary,  especially  in  the  entry  layer.  Accordingly,  the  presumed  diamagnetic  current  (positive 
eastward)  is  opposite  to  the  observed  one.  This  discrepancy  comes  is  due  to  our  assumption  that 
the  extrapolated  relation  Pe  ~  (l/2)Pi  upstream  of  the  SAID  channel  is  still  valid  inside  the 
channel.  However,  while  the  hot  ions  move  farther  inward,  the  hot  electrons  should  accumulate 
about  their  cutoff  boundary.  In  a  steady  state,  this  accumulation  is  balanced  via  precipitation  and 
azimuthal  drift.  The  conjugate  (Polar  and  DMSP)  hot  electron  precipitation  shows  no  noticeable 
increase  near  the  boundary.  Therefore,  the  azimuthal  (diamagnetic)  drift  is  the  prime  contender 
thereby  giving  the  observed  westward  current  in  the  entry  layer  [10]. 
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Figure  5,  Fine  structure  of  the  SAID  channel 

It  is  readily  seen  that  the  suprathermal  ion  population  is  coineident  with  enhaneed  wave  aetivity 
in  the  channel,  and  ~l-eV  ions  extend  also  up-  and  down-stream.  Note  that  in  the  entry  layer  MS 
means  the  ICB  waves.  As  shows  panel  d,  the  inner  SAID  boundary  corresponds  to  the  drop  of  Pi 
and  Ni  to  =0.3  nPa  and  0.4  cm”^,  respectively,  i.e.  by  a  factor  of  4,  which  is  consistent  with  the 
value  ofNmin-  As  the  hot  ion  anisotropy  increases  across  the  channel,  5T  =  0.1  — 0.8,  the  (IB) 
ion  distribution  reduces  to  the  so-called  ion-ring  distribution,  which  is  the  main  driver  of  the 
wave  activity  near  the  center  of  the  channel  and  further  inward.  Note  also  that  dEmis'^ 
downstream  of  the  channel  remains  well  above  the  upstream  level.  This  observation  is  the  key 
for  explaining  the  relation  of  the  outer  RB  boundary  with  the  plasmapause  earlier  discussed  in 
[16]. 

4.3  Physics-based  hybrid  model  with  finite  Larmor  radius  effects 

Strong  UHF/GPS  scintillations  in  the  subauroral  ionosphere  indicate  the  presence  of  density 
irregularities  with  transverse  scale  lengths  <100  m.  As  the  state-of-the-art  fluid  models,  e.g. 
AFRL/PBMOD,  are  incapable  of  predicting  small-scale  structures  ~10s-100s  of  meters,  it  is 
necessary  to  go  beyond  the  fluid  approximation  and  include  ion  kinetic  effects.  A  numerical, 
first-principle  hybrid  model  of  interchange  and  Kelvin-Helmholtz  instabilities  in  the  equatorial 
and  midlatitude  ionosphere  has  been  developed  [17-19].  Nonlinear  system  of  equations 
describing  the  evolution  of  the  interchange  and  Kelvin-Helmholtz  (KH)  instabilities  in  the 
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presence  of  finite  Larmor  radius  effects  in  low-  (ionosphere)  and  high-beta  (plasma  sheet) 
plasmas  with  velocity  shear  has  been  derived.  It  includes  the  combined  effect  of  interchange  and 
KH  instabilities  in  2D  geometry  in  a  plane  perpendicular  to  an  external  magnetic  field  (||Z)  with 
transverse  (IX)  inhomogeneity  in  plasma  density  and  plasma  flow  (||Y).  The  normalized  system 
of  collisionless  fluid  equations  is  as  follows 


dh.5(p  dASo)  „  dSw  T.  ,  d^Sn  T.  „  dSn 

— -  =  -K — -+K  — -+—K - +  K  — 

ot  oy  oy  4  ovqv  oy 

-\-  —  y -{V 5(p,5n]-[5(p,A5(p]-\-V^^  [5n,5(p], 


(1) 


dSn  ddn  f  „  „  , 

—  =  -V,  —  +  {5n,5(p}.  (2) 

dt  dy 

Equation  (1)  describes  the  evolution  of  the  perturbed  vorticityAJ^^  ,  Eq.  (2)  describes  the 
perturbed  density  <?«  ,  and  5(p  is  the  first-order  potential.  Einite  gyro-radius  effects  are  captured 
through  the  terms  that  include  the  ratio  of  the  electron  to  ion  temperatures  (T/T)  •  Equations  (1) 
and  (2)  are  in  dimensionless  form,  with  time  normalized  by  the  ion  cyclotron  frequency  and 
with  spatial  scales  normalized  by  p,  =  E  /D„ ,  where  E  =  is  the  sound  speed  with  M,  the 

ion  mass.  In  Eqs.  (I)  and  (2),  A  =  E  =  dV^jclx,  and  E'  =  d^V^jdx^.  The  Poisson  bracket  operator 
is  defined  as 


dx  dy  dy  dx  ' 


The  y-directed  sheared  velocity  function  selected  for  this  analysis  is 


X-Xg 

L 


(3) 


where  Emax  is  the  peak  velocity  amplitude,  L  determines  the  gradient  of  the  shear,  centered  at  Xg. 

Equations  (1)  and  (2)  are  solved  numerically  on  a  uniform  finite-difference  grid  in  2D  {x,y).  A 
2D  code  previously  developed  for  the  analysis  of  flute-type  unstable  modes  [9]  has  been  adapted 
to  solve  the  nonlinear  equations.  The  spatial  grid  extends  between  l<x,y<  256,  where  x  andy  are 
normalized  spatial  dimensions,  and  dx  =  dy  =  1.  Periodic  boundary  conditions  are  used  in  both 
directions.  The  equations  are  solved  using  a  predictor-corrector  scheme,  with  a  fixed  time  step 
dt  =  0.5  (normalized).  Artificial  damping  and  viscosity  terms  are  used  for  numerical  stability. 

Einear  perturbations  of  the  quantities  Sn  and  d(p  of  the  form 

5p[]  ^/)(x)exp[/(Ev-ft)t)], 
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(4) 


are  considered,  where  <«  =  (»,.  +  iy  is  the  complex  frequency  with  y>0  indicating  linear  instability. 
Substitution  into  Eq.  (2)  as  dropping  the  nonlinear  term  shows  that  SnUQ  in  the  linear  regime. 
From  Eq.  (1),  again  dropping  the  nonlinear  Poisson  bracket  terms,  the  following  second  order 
differential  equation  for  the  linearly  perturbed  potential  S(p  is  obtained; 


d^5(p 

dx^ 


e5(p- 


kv: 


CO-kVr, 


5(p  -  0. 


(5) 


This  equation  is  solved  numerically  on  a  ID  (x)  grid  using  a  standard  shooting  code  with  zero 
derivative  {d5(pldx  =  Q)  boundary  conditions.  Using  L  =  32,  x„  =  129,and  Eomax=h  the  numerical 
solution  to  Eq.  (5)  is  found  and  used  to  obtain  the  growth  rate  as  a  function  k.  This  result  is 
shown  in  Figure  6(a)  as  the  solid  black  curve.  The  red  dot  indicates  the  expected  growth  rate  for 
the  longest  unstable  wavelength  that  fits  with  the  simulation  box;  kL  =  InLjX  □  0.7854.  Figure  6(b) 

plots  the  spatially-averaged  square  of  the  potential  on  the  grid  as  a  function  of  time 

(solid  black  curve).  The  dashed  red  line  is  an  approximate  fit  to  the  black  curve  and  demonstrates 
that  the  calculated  linear  growth  rate  agrees  well  with  that  from  dispersion  analysis  for  the 
Kelvin-Helmholtz  instability. 


Figure  6,  Linear  dispersion  analysis 

Then,  a  2D  numerical  code  was  run  with  these  same  parameters  to  obtain  the  linear  growth  of  the 
Kelvin-Helmholtz  instability  and  nonlinear  saturation  (Figure  7).  The  linear  growth  of  the 
instability  saturates  due  to  the  nonlinear  interaction  of  different  unstable  modes  in  the  system. 
During  this  phase,  several  of  the  nonlinear  terms  in  Eqs.  (1)  and  (2)  provide  coupling  between 
the  potential  and  density.  Eventually,  broad-spatial-spectrum  density  structures  appear. 
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Figure  7,  Time  evolution  of  the  potential 

Figure  8  shows  sample  results  for  a  eode  run  with  =  0.1,  TjT.=l,L  =  8,  and  Xg  =  129.  The 
averaged  squared  potential  as  a  function  of  time  is  shown  in  Figure  8(a).  The  linear  growth  phase 
occurs  in  the  first  -1000  ion  cyclotron  periods  and  is  not  visible  on  this  scale.  After  the  system 
nonlinearly  saturates,  the  potential  structures  undergo  an  oscillatory  phase  followed  by  a  very 
slowly  evolving  static  state  where  the  average  potential  falls  linearly  over  roughly  3  million  ion 
cyclotron  periods. 


Figure  8,  Nonlinear  regime  of  the  instahility 
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After  □  3.5  X 10'’ ,  the  system  begins  to  rapidly  evolve,  with  nonlinear  coupling  driving  density 
structures.  The  (first-order)  potential  and  density  at  Qj  □  3.7  xio'’  on  the  grid  is  shown  in  frames 
(b)  and  (c).  The  density  contours  in  frame  (c)  show  evidence  of  long  and  short  wavelength 
features.  To  illustrate  this,  frame  (d)  plots  the  spectrum  for  the  density  at  this  time  (averaged  in 
they  direction).  This  spectrum  is  substantially  broader  than  the  density  spectrum  at  earlier 
times  in  this  simulation.  This  demonstrates  the  existence  of  nonlinear  cascade  to  the  short 
wavelength  part  of  the  spectrum.  In  frame  (d)  one  also  can  see  that  the  spectrum  is  of  non- 
Kolmogorov  type  with  different  spectral  indexes  for  short  and  long  wavelength  parts  of  the 
spectrum. 


5.  CONCLUSIONS 

Our  study  addresses  basic  problems  in  the  physics  of  the  subauroral  geospace  associated  with  the 
most  intense  midlatitude  space  weather  events  including  UHF/GPS  and  OTHR  disruptions  that 
indicate  the  excitation  of  short  (ion  gyro  radius-scale)  structures.  We  have  explored 
multispacecraft  magnetically-conjugate  observations  of  substorm  SAID  events  near  the  magnetic 
equator  and  in  the  topside  ionosphere.  The  overall  features  negate  the  paradigm  of  voltage  and 
current  generators  and  agree  with  a  short-circuiting  of  substorm-injected  hot  plasma  jets  over  the 
plasmapasphere  and  formation  of  a  turbulent  boundary  layer.  Enhanced  plasma  turbulence  within 
the  SAID  channel  provides  anomalous  circuit  resistivity  and  magnetic  diffusion,  as  in  the  well- 
documented  plasmoid-magnetic  barrier  problem.  The  suprathermal  plasma  population  is 
enhanced  inside  the  channel  and  earthward.  The  short-circuiting  occurs  when  the  cold  plasma 
density  exceeds  a  critical  value  and  the  polarization  field  at  the  front  is  shorted  out.  As  a  result, 
the  hot  electrons  are  arrested,  while  the  hot  ions  yet  move  inward.  This  provides  a  natural 
explanation  of  the  long-known  dispersionless  auroral  boundary.  However,  we  still  know  little 
about  the  underlying  mechanism  of  SAPS,  particularly  SAPS  wave  structures  (SAPSWS),  except 
that  nonlinear  effects  are  significant.  SAPSWS  appear  with  RB  and  RC  precipitation,  radar 
clutter,  and  UHF/GPS  scintillations  after  the  onset  of  substorms.  The  AFRL  hybrid  code  "Flute" 
has  been  modified  into  a  new  hydrodynamic  code  with  the  ion  gyro  kinetics  for  modeling 
interchange  and  Kelvin-Helmholtz  instabilities  in  the  equatorial  and  midlatitude  ionosphere 
creating  short  scale  turbulence  responsible  for  UHF/GPS-F  band  scintillations.  Nonlinear  system 
of  equations  describing  the  evolution  of  the  instabilities  in  the  presence  of  finite  Farmer  radius 
effects  in  low-  (ionosphere)  and  high-beta  (plasma  sheet)  plasmas  with  velocity  shear  has  been 
derived.  A  long-time-scale  calculation  of  2D  nonlinear  equations  shows  the  formation  of 
nonlinear  density  structures  with  a  broad  wave-number  spectrum.  The  spectrum  is  substantially 
broader  than  the  density  spectrum  at  earlier  times  in  this  simulation.  This  demonstrates  the 
existence  of  nonlinear  cascade  to  the  short-wavelength  part  of  the  spectrum  of  non-Kolmogorov 
type  with  different  spectral  indexes  for  short  and  long  wavelengths. 
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